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ABSTRACT 

Context. Tidal interactions between planets and their host stars are not well understood, but may be an important factor in their 
formation, structure and evolution. Previous studies suggest that these tidal interactions may be responsible for discrepancies between 
the ages of exoplanet host stars estimated using stellar models (isochronal ages) and age estimates based on the stars’ rotation periods 
(gyrochronological ages). Recent improvements in our understanding of the rotational evolution of single stars and a substantial 
increase in the number of exoplanet host stars with accurate rotation period measurements make it worthwhile to revisit this question. 
Aims. Our aim is to determine whether the gyrochronological ages for transiting exoplanet host stars with accurate rotation period 
measurements are consistent with their isochronal ages, and whether this is indicative of tidal interaction between the planets and their 
host stars. 

Methods. We have compiled a sample of 28 transiting exoplanet host stars with measured rotation periods, including two stars (HAT- 
P-21 and WASP-5) for which the rotation period based on the light curve modulation is reported here for the first time. We use our 
recently-developed Bayesian Markov chain Monte Carlo method to determine the joint posterior distribution for the mass and age of 
each star in the sample. We extend our Bayesian method to include a calculation of the posterior distribution of the gyrochronological 
age that accounts for the uncertainties in the mass and age, the strong correlation between these values, and the uncertainties in the 
mass-rotation-age calibration. 

Results. The gyrochronological age (r gyro ) is significantly less than the isochronal age for about half of the stars in our sample. Tidal 
interactions between the star and planet are a reasonable explanation for this discrepancy in some cases, but not all. The distribution 
of r gyro values is evenly spread from very young ages up to a maximum value of a few Gyr, i.e., there is no obvious “pile-up” of stars 
at low values or very high of T gyro as might be expected if some evolutionary or selection effect were biasing the age distribution of 
the stars in this sample.. There is no clear correlation between T gyr0 and the strength of the tidal force on the star due to the innermost 
planet. There is clear evidence that the isochronal ages for some K-type stars are too large, and this may also be the case for some 
G-type stars. This may be the result of magnetic inhibition of convection. The densities of HAT-P-11 and WASP-84 are too large 
to be reproduced by any stellar models within the observed constraints on effective temperature and metallicity. These stars may 
have strongly enhanced helium abundances. There is currently no satisfactory explanation for the discrepancy between the young age 
for CoRoT-2 estimated from either gyrochronology or its high lithium abundance, and the extremely old age for its K-type stellar 
companion inferred from its very low X-ray flux. 

Conclusions. There is now strong evidence that the gyrochronological ages of some transiting exoplanet host stars are significantly 
less than their isochronal ages, but it is not always clear that this is good evidence for tidal interactions between the star and the planet. 

Key words, stars: solar-type - planet-star interactions 


1. Introduction 

are bom rotating rapidly and can then lose angular mo¬ 
mentum if they have a magnetised stellar wind. This observa¬ 
tion leads to methods to estimate the age of single late-type stars 
from their rotation period, a technique known as gyrochronol¬ 
ogy (Barnes 2007). It is unclear whether exoplanet host stars can 
be considered as single stars in this context. If significant orbital 
angular momentum is transferred by tides from the orbit of a 
planet to the rotation of the host star (“tidal spin-up”) then the 
star may rotate faster than a genuine single star of the same age, 
i.e., its gyrochronological age (T gyro ) will be an underestimate of 
its true age. There is currently no quantitative global theory to 
calculate the efficiency of tidal spin-up for stars with convective 
envelopes because the dissipation of tidal energy by the exci¬ 
tation and damping of waves in convective atmospheres is not 


fully understood (Zahn 1975, 1977; Goodman & Lackner2009; 
Damiani & Lanza 2014; Ogilvie 2014). Alternatively, the planet 
may disrupt the magentic field geometry of the star and thereby 
reduce the efficiency of angular momentum loss in the magne¬ 
tised stellar wind, although this is only likely to be a significant 
effect for F-type stars (T e ff > 6000 K) with short period planets 
(P < lOd, Lanza 2010). 

If a planet transits its host star then an analysis of the eclipse 
light curve can yield an accurate estimate for the radius of the 
star relative to the semi-major axis of the planet’s orbit, R*/a, 
provided that the eccentricity of the orbit is known. This esti¬ 
mate can be combined with Kepler’s laws to estimate the density 
of the host star (Seager & Mallen-Ornelas 2003). The density 
can be combined with estimates for the effective temperature 
and metallicity of the star to infer a mass and age for the star 
by comparison with stellar models. The age derived by compar- 
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ing the properties of a star to a grid of stellar models is known 
as the isochronal age (ri s0 ). It is not straightforward to calculate 
the statistical error on Ti s0 because it is strongly correlated with 
the estimate of the stellar mass (M*) and the distribution can be 
strongly non-Gaussian, e.g., for stars near the end of their main- 
sequence lifetime the probability distribution for Ti so given the 
observed constraints (the posterior probability distribution) can 
be bi-modal (Maxted et al. 2015). There are also systematic er¬ 
rors in stellar models due to uncertainties in the input physics. 
For example, the efficiency of energy transport by convection 
is usually parameterised by the mixing length parameter »mlt 
( the path length of a convective cell in units of the pressure scale 
height) that is not known a-priori but that can be estimated by 
finding the value of ctmlt for which the models match the ob¬ 
served properties of the Sun. The initial helium abundance for 
the Sun is estimated in a similar way. The radii of some M- and 
K-type dwarf stars are not accurately predicted by these “solar- 
calibrated” stellar models (Hoxie 1973; Popper 1997). This “ra¬ 
dius anomaly” for low mass stars is an active research topic that 
is motivated by a need to better understand exoplanet host stars 
and is driven by advances in simulating convection in low-mass 
stars (Ludwig et al. 2008) and incorporating magnetic fields into 
stellar models (Feiden & Chaboyer 2013). 

Pont (2009) found suggestive empirical evidence for excess 
rotation in stars with close-in giant planets (“hot Jupiters”) based 
on a sample of 28 stars. For 3 of these stars the rotation period 
was measured from the rotational modulation of the light curve 
by star spots. For the other 25 stars the rotation period was esti¬ 
mated from the projected equatorial rotation velocity measured 
from rotation line broadening (V rot sin/*). For 12 of the stars, 
V ro t sin ;* was below the instrumental resolution. It is difficult to 
make accurate estimates of V rot sin /* in these cases because it 
requires the deconvolution of the instrumental profile from the 
observed spectral line profiles and the results are sensitive to de¬ 
tails of the stellar atmosphere models used to calculate the intrin¬ 
sic stellar line profiles. Brown (2014) used a maximum likehood 
method to estimate the ages of 68 transiting exoplanet host stars 
using 5 different grids of stellar models and also applied four dif¬ 
ferent methods to calculate the age using gyrochronology. The 
rotation periods for 8 of the stars in that study were based on ro¬ 
tational modulation of the light curve, the remainder were esti¬ 
mated from y rot sin /*. Estimates of the rotation period based on 
V ro t sin /* are generally much less precise than those measured 
directly from the lightcurve and can be affected by systematic 
errors when the line broadening due to rotation is small com¬ 
pared to other poorly-understood line-broadening effects such as 
macro-turbulence. Additional uncertainty is introduced by the 
unknown inclination of the star’s rotation axis, (*. Brown found 
a “slight tendency for isochrones to produce older age estimates” 
but that the “evidence for any bias on a sample-wide level is in¬ 
conclusive.” 

Planet host stars are an interesting test case for our under¬ 
standing of the tidal interactions because tidal spin-up will in¬ 
crease the star’s magnetic activity, which may increase the loss 
of rotational angular momentum through a magnetised stellar 
wind. The balance of these two effects may lead to a quasi¬ 
stationary state in which the planet’s orbit is stable in the long 
term (Damiani & Lanza 2015). However, if tidal interactions 
between stars and planets are strong enough to efficiently de¬ 
stroy massive, short-period planets then we might expect to only 
find hot Jupiters around young stars. To test this hypothesis it 
would be useful in these cases to be able to estimate an accurate 
isochronal age, because it is likely that the star will be spun-up 
during the destruction of the planet and so its gyrochronological 


age will not be an accurate estimate of its true age. Comparing 
the isochronal and gyrochronological ages is not a straightfor¬ 
ward problem because the uncertainties on these two age esti¬ 
mates are strongly non-Gaussian and r gyro depends on the as¬ 
sumed stellar mass, so Ti S0 and r gyro are also correlated. 

In Maxted et al. (2015) we developed a Markov chain Monte 
Carlo (MCMC) method that enables us to estimate the joint 
posterior probability distribution for Ti s0 and M* based on the 
observed density, effective temperature (T e ff) and metallicity 
([Fe/H]) of a star. In this paper we select a sample of 28 tran¬ 
siting exoplanet host stars with accurately measured rotation pe¬ 
riods (Section 2). For two of the stars in this sample the rotation 
periods are new estimates based on the rotational modulation of 
the light curves. We then describe how we have extended the 
method of Maxted et al. (2015) to include a calculation of the 
posterior probability distribution for r gyro and present the results 
of applying this technique to our sample of 28 planet host stars 
(Section 3). The implication of these results for our understand¬ 
ing of tidal spin-up of exoplanet host stars and the reliability of 
solar-calibrated stellar models is discussed in Section 4 and our 
conclusions are given in Section 5. 


2. Sample definition and data selection 

We have compiled a list of 28 transiting exoplanet host stars for 
which the rotation period has been measured directly from the 
rotational modulation of the light curve due to star spots or, in 
the case of 55 Cnc, spectroscopic variability in the Ca n H and K 
emission line fluxes. For all these stars the existence of a transit¬ 
ing planet has been confirmed by multiple radial velocity mea¬ 
surements based on high-resolution spectroscopy or dynamical 
analysis of transit timing variations in multi-planet systems (e.g., 
Kepler-30). The observed properties of these stars are given in 
Table 1. For each star we list the rotation period, P rot , the orbital 
period, P or b, the effective temperature, T e ff, the logarithmic sur¬ 
face iron abundance relative to the Sun, [Fe/H], the mean stellar 
density, />*, the flux from the star at the top of the Earth’s atmo¬ 
sphere, /®, and the logarithm of the stellar luminosity, log(L*). 

We have used only one or two sources of input data for this 
analysis rather than attempting to reconcile or take the average 
of multiple studies for each system. We have excluded stars with 
masses that are too low to be covered by our grid of stellar mod¬ 
els (< 0.6 M q ) or too high for gyrochronology to be applicable 
(> 1.25M 0 ). We have also excluded stars where the density es¬ 
timate is based on a Kepler long-cadance (LC) light curve only, 
because these density estimates appear to be biased compared 
to the density estimated using asteroseismology (Huber et al. 
2013). 

2.1. Luminosity measurements 

We can include the observed luminosity of the star (L*) as an 
additional constraint in the analysis to derive the mass and age 
of the star. For stars that have a trigonometrical parallax in van 
Leeuwen (2007) with precision cr n ln S 0.1 we calculate L* us¬ 
ing the flux at the top of the Earth’s atmosphere (/®) estimated 
by integrating a synthetic stellar spectrum fit by least-squares to 
the observed fluxes of the star. Optical photometry is obtained 
from the Naval Observatory Merged Astrometric Dataset (NO¬ 
MAD) catalogue 1 (Zacharias et al. 2004), the Tycho-2 catalogue 
(H0g et al. 2000), The Amateur Sky Survey (TASS, Droege et al. 


1 http://www.nofs.navy.mi1/data/fchpix 
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Table 2. Periodogram analysis of the WASP light curves for WASP-5. 
Observing dates are JD-2450000, N is the number of observations used 
in the analysis, A is the semi-amplitude of the sine wave fit by least- 
squares at the period P found in the periodogram with false-alarm prob¬ 
ability FAP. 


Camera 

Dates 

N 

B[d] 

A [mag] 

FAP 

227 

3870-4054 

4429 

8.18 

0.003 

0.027 

228 

3870-4020 

3918 

16.35 

0.003 

0.094 

228 

4284-4433 

4673 

8.52 

0.002 

0.052 

225 

5352-5527 

7720 

15.04 

0.006 

< 0.001 

225 

5716-5897 

4866 

16.25 

0.002 

0.001 

228 

5731-5885 

3954 

18.64 

0.003 

0.213 


1997) and the Carlsberg Meridian Catalog 14 (Copenhagen Uni¬ 
versity et al. 2006). Near-infrared photometry was obtained from 
the Two Micron All Sky Survey (2MASS) 2 and Deep Near In¬ 
frared Survey of the Southern Sky (DENIS) 3 catalogues (The 
DENIS Consortium 2005; Skrutskie et al. 2006). The synthetic 
stellar spectra used for the numerical integration of the fluxes 
are from Kurucz (1993). Reddening can be neglected for these 
nearby stars given the accuracy of the measured fluxes and paral¬ 
laxes. Standard errors are estimated using a simple Monte Carlo 
method in which we generate 65,536 pairs of n and /® values 
from Gaussian distributions and then find the 68.3% confidence 
interval of the resulting log(L/L®) values. The 2MASS photom¬ 
etry for 55 Cnc is not reliable so for this star we have used the 
value of /® from Boyajian et al. (2013). 

2.2. Rotation periods 

The rotation periods for the stars in Table 1 are taken from the 
published sources noted except for HAT-P-21 and WASP-5, for 
which we have measured the rotation periods using observa¬ 
tions obtained by the WASP project (Pollacco et al. 2006) and 
the method described by Maxted et al. (2011). For HAT-P-21 
we used 4185 observations obtained between 2007 Jan 2 and 
2007 May 15. These data show a very clear sinusoidal signal 
with a period of 15.9 d with an amplitude of 0.012 magnitudes. 
This period agrees very well with the estimate P mt = 15.7 + 2.2 d 
that is expected based on the projected equatorial rotation veloc¬ 
ity (V r otsin;* = 3.5 ± 0.5km s^ 1 , Bakos et al. 2011) and the 
radius of the star, assuming that the inclination of the star’s rota¬ 
tion axis is i ~ 90°. The WASP data for HAT-P-21 are shown as 
a function of rotation phase in Fig. 1. The standard error on this 
value has been estimated from the full-width at half-maximum 
of the peak in the periodogram. 

For WASP-5 we analysed 6 sets of data obtained with 3 dif¬ 
ferent cameras in 4 different observing seasons, as detailed in 
Table 2. In all 6 data sets we detect a periodic signal consistent 
with a rotation period P rot ~ 16 d if we allow for the possibility 
that the distribution of star spots can produce a signal at the first 
harmonic of the rotation period (i.e., at P mt /2). The value of P I0t 
in Table 1 is the mean of the observed values calculated with 
this assumption and the error quoted in the standard error of the 
mean. 

Mohler-Fischer et al. (2013) report a factor of 2 ambiguity 
in the rotational period of HATS-2. We have used P rot = 24.98 d 
for our analysis since this gives a value of r gyro that is more con¬ 
sistent with the value of r; so . 

2 http://www.ipac.caltech.edu/2mass 

3 http://cdsweb.u-strasbg.fr/denis.html 



Phase 


Fig. 1. WASP photometry for HAT-P-21 plotted as a function of rotation 
phase assuming P rot = 15.9 d. 


3. Analysis and Results 

We used version 1.1 of the program bagemass 4 (Maxted et al. 
2015) to calculate the joint posterior distribution for the mass 
and age of each star based on the observed values of T e ff, [Fe/H], 
the mean stellar density p* and, if available, L*. The stellar mod¬ 
els used for our analysis were produced with the garstec stellar 
evolution code (Weiss & Schlattl 2008). The initial composition 
of the models is computed assuming a cosmic helium-to-metal 
enrichment AY/AZ = (7® - Ibbn)/Z®, where Tbbn = 0.2485 
is the primordial helium abundance due to big-bang nucleosyn¬ 
thesis (Steigman 2010), Z® = 0.01826 is the solar metal abun¬ 
dance, and the initial solar helium abundance is 7® = 0.26646 so 
AY/AZ = 0.984. These models do not include rotation, but the 
rotation rates of the stars in our sample are low (<2% of breakup 
velocity in all but two cases) and so it is not expected that rota¬ 
tion would play a direct role in the structure of these stars. We do 
not rule out other, indirect effects, e.g. low convective efficiency 
associated to large magnetic fields, but these have to be modelled 
in an ad-hoc manner. This is discussed later in this work, in rela¬ 
tion to specific stars in the sample. The methods used to calculate 
and interpolate the stellar model grid are described in Serenelli 
et al. (2013) and Maxted et al. (2015). We set lower limits of 
80 K on the standard error for T e n and 0.07 dex for the standard 
error on [Fe/H] (Bruntt et al. 2010) and assume flat prior dis¬ 
tributions for the stellar mass and age. The results are given in 
Table 3, where the maximum-likelihood (best-fit) values of the 
stellar mass and age are denoted Mb and r lso |-,, respectively. The 
mean and standard deviation of the posterior distributions for the 
mass and age are listed under (M*) and (ri so ). Also listed in this 
table are our estimates of the systematic errors in these values 
due to an assumed error of 0.2 for (Tmlt and an assumed error 
of 0.02 for the initial helium abundance, 7. The change in the 
estimated mass and age due to increasing 7 by its assumed error 
are given by the quantities cr M Y and cr T j , respectively. Similarly, 
c tmm and cr T a quantify the change in the estimated mass and age 
due to the error in crMur- We show the best-fit values of M* and 
Tiso and the effects of changing 7 and o-vna by their assumed 
uncertainties for all the stars in the sample in Fig. 2. 


4 http://sourceforge.net/proj ects/bagemass 
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Table 1. Observed properties of stars in our sample. 


Star 

Prot 

Porb 

T e ff 

[Fe/H] 

P* 

/© 

log(L*) 

Ref. 


Id] 

Id] 

[K] 


[Po] 

[pWm' 2 ] 

[L 0 ] 


55 Cnc 

39.00 ± 9.00 

0.74 

5234 ± 30 

+0.31+0.04 

1 -^o4_0.036 

120.4+1.0 

-0.244 ± 0.009 

1,2, 3,4 

CoRoT-2 

4.52 ± 0.02 

1.74 

5598 ± 50 

+0.04 ± 0.05 

1.362 + 0.064 



5,6 

CoRoT-4 

8.87+1.12 

9.20 

6190 + 60 

+0.05 ± 0.07 

0 790 +0064 
U - ,7U -0.161 



7,8 

CoRoT-6 

6.40 ± 0.50 

8.89 

6090 ± 70 

-0.20 + 0.10 

0.929 ± 0.064 



9,8 

CoRoT-7 

23.64 + 3.62 

0.85 

5313 + 73 

+0.03 ± 0.07 

1.671 +0.073 



10, 11, 12 

CoRoT-13 

13 00 +5 °° 

1J.UU_ 300 

4.04 

5945 ± 90 

+0.01 ± 0.07 

0.526 ± 0.072 



13,8 

CoRoT-18 

5.40 + 0.40 

1.90 

5440 ±100 

-0.10 + 0.10 

1.090 + 0.160 



14,6 

HAT-P-11 

30.50^0 

4.89 

4780 ± 50 

+0.31+0.05 

2.415 ± 0.097 

5.9 ±0.3 

-0.590 ± 0.035 

15,8 

HAT-P-21 

15.90 + 0.80 

4.12 

5634 ± 67 

+0.04 ± 0.08 

0.700 + 0.150 



16, 12 

HATS-2 

24.98 ± 0.04 

1.35 

5227 ± 95 

+0.15+0.05 

1.220 + 0.060 



17, 17 

HD 189733 

11.95 + 0.01 

2.22 

5050 ± 50 

-0.03 ± 0.05 

1.980 + 0.170 

27.5 ± 1.4 

-0.489 ± 0.024 

18, 19 

HD 209458 

10.65 ± 0.75 

3.52 

6117 + 50 

+0.02 ± 0.05 

0.733 ± 0.008 

23.1 ± 1.2 

0.248 ± 0.041 

20, 19 

Kepler-17 

12.10+1.56 

1.49 

5781 +85 

+0.26 + 0.10 

1 121 +0 - 015 
A,1Z ' A -0.034 



21,6 

Kepler-30 

16.00 ± 0.40 

29.33 

5498 ± 54 

+0.18+0.27 

1.420 + 0.070 



22, 23, 22 

Kepler-63 

5.40 + 0.01 

9.43 

5576 ± 50 

+0.05 ± 0.08 

1 "MS +0.089 
-0.083 



24, 24 

Qatar-2 

11.40 + 0.50 

1.34 

4645 ± 50 

-0.02 ± 0.08 

1.591+0.016 



25,25 

WASP-4 

22.20 + 3.30 

1.34 

5540 ± 55 

-0.03 ± 0.09 

1.230 + 0.022 



26,6 

WASP-5 

16.20 ± 0.40 

1.63 

5770 ± 65 

+0.09 ± 0.09 

0.801 ± 0.080 



6 

WASP-10 

11.91+0.05 

3.09 

4675 ± 100 

+0.03 ± 0.20 

2 7S9 +0 053 



27, 28, 29 

WASP-19 

11.76 + 0.09 

0.79 

5460 ± 90 

+0.14 + 0.11 

0.885 ± 0.006 



30,31 

WASP-41 

18.41+0.05 

3.05 

5450 ± 150 

-0.08 ± 0.09 

1.270 + 0.140 



32, 32 

WASP-46 

16.00+1.00 

1.43 

5600 ±150 

-0.37 + 0.13 

1.240 + 0.100 



33,33 

WASP-50 

16.30 + 0.50 

1.96 

5400 ±100 

-0.12 + 0.08 

1.376 + 0.032 



34, 35, 34 

WASP-69 

23.07 + 0.16 

3.87 

4700 ± 50 

+0.15+0.08 

1.540 + 0.130 



36, 36 

WASP-77 

15.40 + 0.40 

1.36 

5500 ± 80 

+0.00 + 0.11 

1 1 ^7+0.016 

1 * 1 7 -0.020 



37, 37 

WASP-84 

14.36 + 0.35 

8.52 

5300 ±100 

+0.00 + 0.10 

2.015 ± 0.070 



36, 36 

WASP-85 

14.64+1.47 

2.66 

5685 ± 65 

+0.08+0.10 

1.280 + 0.010 



38,38 

WASP-89 

20.20 ± 0.40 

3.36 

4960 ±100 

+0.15+0.14 

i oc 7+0.069 
-0.075 



39, 39 


References. (1) Henry et al. (2000); (2) Dragomir et al. (2014); (3) Valenti & Fischer (2005); (4) Boyajian et al. (2013); (5) Lanza et al. (2009); 
(6) Southworth (2012); (7) Aigrain et al. (2008); (8) Southworth (2011); (9) Fridlund et al. (2010); (10) Lanza et al. (2010); (11) Barros et al. 
(2014); (12) Torres et al. (2012); (13) Cabrera et al. (2010); (14) Hebrard et al. (2011); (15) Sanchis-Ojeda & Winn (2011); (16) Bakos et al. 
(2011); (17) Mohler-Fischer et al. (2013); (18) Henry & Winn (2008); (19) Southworth (2010); (20) Silva-Valio (2008); (21) Beky et al. (2014); 
(22) Sanchis-Ojeda et al. (2012); (23) Fabrycky et al. (2012); (24) Sanchis-Ojeda et al. (2013); (25) Mancini et al. (2014); (26) Sanchis-Ojeda 
et al. (2011); (27) Smith et al. (2009); (28) Christian et al. (2009); (29) Barros et al. (2013); (30) Tregloan-Reed et al. (2013); (31) Mancini et al. 
(2013); (32) Maxted et al. (2011); (33) Anderson et al. (2012); (34) Gillon et al. (2011); (35) Tregloan-Reed & Southworth (2013); (36) Anderson 
et al. (2014); (37) Maxted et al. (2013); (38) Brown et al. (2014); (39) Hellier et al. (2014). 


3.1. Gyrochronological ages 

Our Bayesian MCMC method produces a large set of points in 
the mass-age parameter space (“Markov chain”) that has the 
same distribution as the posterior probability distribution for 
these parameters. We used equation (32) from Barnes (2010) to 
calculate T gyro from the rotation period and the mass of the star 
for every point in the Markov chain for each star. The convec¬ 
tive turn-over time scale that encapsulates the mass dependence 
of this mass-age-rotation relation was interpolated from Table 
1 of Barnes & Kim (2010). This table suggests that stars with 
masses M* 1 ,25M e do not have convective envelopes so we 
have restricted our analysis to stars with P(M* < 1.25M 0 ) > 0.9 


according to our MCMC results. We also checked that none of 
the stars are close to the end of the main sequence, where rapid 
changes in radius and internal structure are likely to make gy¬ 
rochronological ages unreliable. The calculation of r gyr0 requires 
an estimate for the value of the parameter Pq. We account for the 
uncertainty in this value by randomly generating a value of Pq 
uniformly distributed in the range 0.12 d to 3.4 d for each point 
in the Markov chain. In principle, there is some additional uncer¬ 
tainty in T gyro due to the stars’ surface differential rotation com¬ 
bined with the variation in the latitude of the active regions that 
produce the star spot modulation (Epstein & Pinsonneault 2014). 
We assume that this uncertainty affects the calibration sample 
used by Barnes (2010) to the same extent that it affects the stars 
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Fig. 2. Change in the best-fitting masses and ages of transiting exoplanet host stars due to a change in the assumed helium abundance or mixing 
length parameter. Dots show the best-fitting mass and age for the default values of Y and o:mlt and lines show the change in mass and age due to 
an increase in helium abundance AY = +0.02 (left panel) or a change in mixing length parameter A<r MLT = -0.2 (right panel). Horizontal lines 
indicate the age of the Galactic disc (dashed), the age of the Universe (dotted) and the largest age in our grid of stellar models (solid). The curved 
dotted line shows the terminal age main sequence (TAMS) for stars with solar composition. 


in our sample so that this uncertainty is already accounted for 
by randomly perturbing the parameter Pq. The uncertainties in 
the observed values of P TOt are accounted for in a similar way, 
but using Gaussian random distributions for the standard errors 
shown in Table 1. Fig. 3 shows the joint posterior distributions 
for (M*, ri S o) and (t; so , Tg yro ) calculated from the Markov chain 
using this method for CoRoT-13. 

For each star we compare the values of r gyro calculated using 
the mass value at every point in the Markov chain to the cor¬ 
responding Ti S o values for the same point in the Markov chain. 
To quantify the difference between t; so and r gyro we calculate 
the fraction of points p T in the chain for which r gyro > Ti S0 , i.e., 
p T — P(r gyro > Tj s0 ) is the probability that the gyrochronogical 
age is greater than the isochronal age. The results of this com¬ 
parison are given in Table 4 and are shown in Fig. 4. Also given 
in Table 4 is Ttidai> which is a very approximate estimate of the 
time scale for tidal interactions between the star and the planet 
(Albrecht et al. 2012). It must be emphasized that the actual time 
scale for tidal spin-up in these systems is uncertain by a few or¬ 
der of magnitude (Ogilvie 2014). For multi-planet systems, r t id a i 
applies to the innermost planet. The values of r gyro as a function 
of t tidal are shown in Fig. 5. 

4. Discussion 

It is clear that the gyrochronological age is significantly less than 
the isochronal age for about half of the stars in this sample. This 
discrepancy is apparent from Fig. 4, but the errors on these val¬ 
ues are correlated and non-Gaussian so to accurately quantify 
this discrepancy we need to use the values of p T in Table 4. There 
is no obvious relation between the gyrochronological age of the 
stars and the estimated time scale for tidal interactions between 
the star and the planet. This may not be surprising given that 
t ddai is uncertain by a few orders of magnitude (Ogilvie 2014). It 
is also apparent that there is an even spread of r gyr0 values from 
very young ages up to a maximum value of a few Gyr, i.e., there 
is no obvious “pile-up” of stars with very low or very high r gyro 
values as might be expected if some evolutionary or selection ef¬ 
fect were biasing the age distribution of the stars in this sample. 

The discrepancy between isochronal and gyrochronological 
ages for some planet host stars has been noted for more lim¬ 
ited samples of planet host stars before, and has been cited as 



Tso [Gyr] 


Fig. 4. Comparison of gyrochronological ages (r gyro ) to isochronal ages 
(r is0 ) for planet host stars with measured rotation periods. Points with 
error bars indicate the mean and standard deviation of the posterior age 
distribution. The straight line is the relation r gyro = r iso . 


evidence for tidal spin-up of the stars by their planetary com¬ 
panions, or for the disruption of the normal spin-down process 
for these stars (Pont 2009; Lanza 2010; Poppenhaeger & Wolk 
2014). However, there are other possibile explanations in some 
cases. It may be that the isochronal ages for some stars are not 
accurate because of missing physics in the stellar models. It is 
also possible for the distribution of active regions on a star to 
produce a modulation of the light curve with a period of half 
the rotation period. This is not a viable explanation for the low 
values of p T in Table 4 in general, but may be an issue for a 
few stars with limited photometric data. We have been careful to 
select stars for our sample for which the observational data are 
robust, so the general appearance of distribution in Fig. 4 can¬ 
not be ascribed to observational errors, although there may be 
issues with the observed data for a few stars. The limitations of 
the observed data, the stellar models and the estimated value of 
t udai affect different stars in this sample in different ways, so to 
understand the implication of these results we need to look at the 
results on a case-by-case basis. 
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Table 3. Bayesian mass and age estimates for the host stars of transiting extrasolar planets using garstec stellar models assuming ar MLT = 1-78. 
Columns 2, 3 and 4 give the maximum-likelihood estimates of the age, mass, and initial metallicity, respectively. Column 5 is the chi-squared 
statistic of the fit for the parameter values in columns 2, 3, and 4. Columns 6 and 7 give the mean and standard deviation of their posterior 
distributions. Column 8 (pms) is the probability that the star is still on the main sequence. The systematic errors on the mass and age due to 
uncertainties in the mixing length and helium abundance are given in columns 9 to 12. 


Star 

Tiso.b [Gyr] 

M h [M e \ 

[Fe/H] ib 

X 1 

<Ti SO > [Gyr] 

<M*> [M 0 ] 

0"r,r 

0- T ,a 

C M.Y 

m M.a 

55 Cnc 

10.9 

0.91 

+0.378 

0.63 

10.91 ± 1.62 

0.913 + 0.020 

0.48 

3.19 

-0.033 

-0.027 

CoRoT-2 

1.8 

0.97 

+0.057 

0.02 

2.66 + 1.62 

0.962 ± 0.034 

1.24 

1.31 

-0.047 

-0.019 

CoRoT-4 

1.2 

1.17 

+0.076 

0.01 

2.09 + 1.07 

1.174 + 0.044 

0.42 

0.74 

-0.046 

-0.026 

CoRoT-6 

3.2 

1.02 

-0.160 

0.01 

3.40+1.49 

1.023 + 0.048 

0.39 

1.68 

-0.037 

-0.035 

CoRoT-7 

1.3 

0.90 

+0.046 

0.02 

2.92 + 1.87 

0.884 ± 0.029 

1.87 

0.91 

-0.045 

-0.004 

CoRoT-13 

6.1 

1.08 

+0.083 

0.02 

5.99 + 1.40 

1.089 + 0.053 

0.30 

1.54 

-0.046 

-0.034 

CoRoT-18 

11.9 

0.85 

-0.011 

0.01 

10.69 + 3.27 

0.868 ± 0.043 

1.37 

3.08 

-0.046 

-0.026 

HAT-P-11* 

0.0 

0.82 

+0.238 

9.26 

0.72 + 0.83 

0.813 + 0.016 

0.00 

-0.00 

-0.023 

0.001 

HAT-P-21 

9.9 

0.96 

+0.121 

0.01 

9.52 + 2.26 

0.971 +0.045 

0.66 

2.61 

-0.045 

-0.035 

HATS-2 

10.0 

0.89 

+0.214 

0.01 

9.70 ± 2.77 

0.892 ± 0.037 

0.99 

3.52 

-0.037 

-0.032 

HD 189733 

1.6 

0.82 

-0.014 

0.02 

4.75 + 3.15 

0.805 ± 0.023 

2.45 

1.09 

-0.036 

-0.007 

HD 209458 

2.4 

1.14 

+0.065 

0.20 

2.42 ± 0.80 

1.143 + 0.038 

0.30 

1.22 

-0.043 

-0.029 

Kepler-17* 

0.2 

1.11 

+0.252 

0.04 

1.48 ± 1.07 

1.075 + 0.034 

0.53 

0.10 

-0.040 

0.033 

Kepler-30 

0.2 

1.01 

+0.192 

0.01 

4.38 + 3.24 

0.919 + 0.065 

1.23 

0.16 

-0.049 

0.018 

Kepler-63 

2.2 

0.97 

+0.074 

0.00 

3.16 + 1.88 

0.958 ± 0.037 

1.08 

1.58 

-0.043 

-0.023 

Qatar-2* 

17.5 

0.75 

+0.149 

3.56 

15.72 + 1.36 

0.767 + 0.013 

0.00 

0.02 

-0.018 

0.013 

WASP-4 

6.2 

0.92 

+0.019 

0.04 

6.27+2.34 

0.919 + 0.045 

0.95 

2.85 

-0.041 

-0.031 

WASP-5 

5.8 

1.03 

+0.150 

0.00 

5.84+1.86 

1.035 + 0.048 

0.62 

2.25 

-0.047 

-0.038 

WASP-10 

0.3 

0.76 

+0.033 

0.03 

6.00 + 4.12 

0.709 + 0.031 

2.56 

0.18 

-0.039 

0.009 

WASP-19 

10.0 

0.94 

+0.222 

0.00 

9.95 ± 2.49 

0.949 ± 0.048 

0.88 

3.20 

-0.044 

-0.039 

WASP-41 

8.4 

0.87 

-0.022 

0.01 

8.25 + 3.59 

0.877 ± 0.049 

1.21 

3.40 

-0.039 

-0.034 

WASP-46 

12.1 

0.80 

-0.294 

0.01 

10.03 + 3.51 

0.833 + 0.051 

0.50 

3.59 

-0.030 

-0.034 

WASP-50 

8.6 

0.85 

-0.054 

0.02 

8.57 + 2.86 

0.851 +0.041 

1.15 

3.33 

-0.038 

-0.030 

WASP-69* 

17.5 

0.76 

+0.251 

0.13 

15.20+1.55 

0.776 ± 0.019 

-0.04 

2.94 

-0.025 

-0.015 

WASP-77 

7.8 

0.91 

+0.055 

0.04 

7.57 + 2.53 

0.918 + 0.049 

0.46 

3.10 

-0.031 

-0.034 

WASP-84* 

0.1 

0.86 

-0.072 

0.73 

1.89 + 1.61 

0.826 ± 0.025 

-0.09 

0.08 

-0.010 

0.013 

WASP-85 

1.0 

1.02 

+0.104 

0.02 

2.10+1.38 

0.994 ± 0.036 

1.12 

0.71 

-0.049 

0.003 

WASP-89 

14.9 

0.81 

+0.236 

0.01 

12.07 + 3.11 

0.844 ± 0.037 

1.00 

4.54 

-0.033 

-0.030 


*Best-fit is for age near the edge of the stellar model grid - cr T y and cr T Q , may not be reliable. 



Fig. 3. Left panel: Joint posterior distribution for the mass and age of CoRoT-13 estimated by our isochrone fitting technique. Right panel: Joint 
posterior distribution for the age of CoRoT-13 estimated by our isochrone fitting technique and using gyrochronology. For clarity, only 10% of the 
points from the Markov chain are plotted. 


4.1. HAT-P-11 and WASP-84 - helium rich stars? 

There is no satisfactory fit to the observed properties of HAT- 
P-11 at any age or mass for our grid of standard stellar models, 
i.e. models using a mixing length calibrated on the solar model 
(Q'mlt = 1.78) and AT/AZ = 0.984, consistent with the pri¬ 


mordial helium abundance and the assumed initial solar com¬ 
position. In Table 5 we summarize the published density esti¬ 
mates for HAT-P-11, including four studies based on the very 
high quality Kepler short-cadence (SC) data for this star. The 
Kepler data clearly show distortions to the light curve due to the 
planet crossing dark spots on the face of the host star. Differ- 


Article number, page 6 of 1 lpage. 11 





P. F. L. Maxted et al.: Ages of transiting exoplanet host stars 


Table 4. Stellar ages measured using our stellar models, r iso , and using gyrochronology, T gyro , and the probability p T = P(T gyro > r iso ). A very 
approximate estimate for the time scale for tidal interaction between the star and the inner-most planet is given under log(T t jd a i). Notes include 
constraints on the age of the star and any companion stars front the X-ray luminosity (r x ), and age constraints from the star’s surface lithium (Li) 
abundance Tu 


Star 

^"iso 

Tgyro 

Pt 

log(T tida i) Notes 


[Gyr] 

[Gyr] 


[yr] 


55 Cnc 

10.91 ± 1.62 

8.10 + 3.54 

0.23 

12.8 

CoRoT-2 

2.66 ± 1.62 

0.17 + 0.06 

0.02 

10.3 

CoRoT-4 

2.10 ± 1.06 

1.81 ± 1.17 

0.39 

14.3 

CoRoT-6 

3.40+1.49 

0.35 + 0.10 

0.01 

13.0 

CoRoT-7 

2.92+1.87 

2.80 + 0.82 

0.52 

13.1 

CoRoT-13 

5.99+1.40 

2.34+1.92 

0.08 

11.9 

CoRoT-18 

10.69 + 3.28 

0.22 ± 0.08 

0.00 

10.2 

HAT-P-11 

0.72 + 0.83 

3.89 + 0.89 

0.98 

15.7 

HAT-P-21 

9.52 + 2.26 

1.64 + 0.29 

0.00 

11.1 

HATS-2 

9.70 + 2.77 

3.10 + 0.30 

0.01 

10.5 

HD 189733 

4.75 + 3.15 

0.71+0.09 

0.07 

11.9 

HD 209458 

2.42 + 0.79 

1.83 + 0.85 

0.29 

12.5 

Kepler-17 

1.48 + 1.07 

1.43 + 0.43 

0.53 

10.3 

Kepler-30 

4.38 + 3.24 

1.47 + 0.24 

0.22 

19.0 

Kepler-63 

3.16+1.88 

0.23 ± 0.06 

0.02 

15.8 

Qatar-2 

15.72+1.36 

0.64 + 0.10 

0.00 

10.2 

WASP-4 

6.27 + 2.34 

2.72 + 0.83 

0.09 

10.6 

WASP-5 

5.84+1.86 

2.13 + 0.52 

0.05 

10.5 

WASP-10 

6.00 + 4.12 

0.66 + 0.10 

0.06 

11.6 

WASP-19 

9.95 + 2.49 

0.89 + 0.12 

0.00 

9.5 

WASP-41 

8.25 + 3.59 

1.71+0.21 

0.03 

12.3 

WASP-46 

10.03 + 3.51 

1.23 + 0.20 

0.01 

10.3 

WASP-50 

8.57 + 2.86 

1.30 + 0.15 

0.00 

11.2 

WASP-69 

15.20+1.55 

2.09 + 0.12 

0.00 

13.9 

WASP-77 

7.57 + 2.53 

1.35 + 0.18 

0.00 

10.3 

WASP-84 

1.89+1.61 

0.99 + 0.10 

0.36 

14.7 

WASP-85 

2.09+1.37 

1.50 + 0.33 

0.41 

12.0 

WASP-89 

12.07 + 3.11 

1.88 + 0.18 

0.00 

10.9 


Companion star Tx > 5 Gyr. All age estimates consistent. 
tu « 0.1 Gyr, t x ~ 0.25 Gyr, but companion rx > 5 Gyr. 
Ti sc , consistent with r gyro but large relative error on both. 
Fast rotator. 

Ti S o consistent with r gyro but large relative error on both. 
Tidal spin-up? 

Tidal spin-up? 

Poor isochrone fit - helium-rich? 

r x = 1 -2 Gyr, companion r x > 5 Gyr => tidal spin-up. 

Tidal spin-up? 

t x — 1 -2 Gyr, companion r x > 5 Gyr => tidal spin-up. 
Ti S o consistent with r gyro but large relative error on both. 
Ti S o consistent with r gyro but large relative error on both. 
Ti S o consistent with r gyro but large relative error on both. 
Fast rotator. 

Inflated K-dwarf. 

Tidal spin-up? 

Tidal spin-up? 

Tidal spin-up? 

Tidal spin-up? 

Magnetically active G-type star - Ti so unreliable? 

Tidal spin-up? 

Tidal spin-up? 

Inflated K-dwarf. 

Companion V rot sin /* => age ~ 0.4 Gyr. 

Magnetically active G-type star - Ti so unreliable? 

Poor isochrone fit - helium-rich? 

Ti S o consistent with T gyro but large relative error on t; so . 
Tidal spin-up? 


5 



8 10 12 14 16 18 20 

| og(T tid0l /yr) 

Fig. 5. Gyrochronological ages as a function of log(T t M a i), a very ap¬ 
proximate estimate of the time scale for tidal interaction between the 
star and the innermost planetary companion. Note that the time scale 
for tidal interactions is uncertain by a few orders of magnitude. Sys¬ 
tems with significant differences between Ti S0 and T gyro (p < 0.05) are 
plotted with filled symbols. 55 Cnc (T iso « 9 Gyr, log(r tidal /y) « 13) is 
not shown here. 


ent investigators have accounted for these spot-crossing events 
in different ways, but the stellar densities derived are all consis¬ 
tent with the value we have used in our analysis. There are also 
three independent analyses of the spectrum of this star that have 
been used to derive the T e ff and [Fe/H]. Again, all three values 
are consistent with each other, including the value we have used 
for our analysis. 

In principle, the Kepler short-cadence photometry for HAT- 
P-l 1 can be used to estimate the density of this star using aster- 
oseismology. Christensen-Dalsgaard et al. (2010) analysed the 
Kepler data for HAT-P-11 from the commissioning period and 
the first month of regular observations and claimed definite evi¬ 
dence for solar-like oscillations, yielding a preliminary estimate 
of its mean density. However, an analysis of the complete Kepler 
data set for this star shows no convincing evidence for solar-like 
oscillations because HAT-P-11 has a much higher level of pho¬ 
tometric noise than a typical star of the same magnitude (Davies, 
priv. comm.). 

In general, the analysis of the light curve for a transiting ex¬ 
trasolar planet makes the assumption that, apart from the effects 
of limb-darkening, the mean surface brightness of a star is the 
same as the surface brightness in the regions obscured by the 
planet. This would not be the case if, for example, the planet 
transits a chord near and approximately parallel to the stellar 
equator of a star with dark spots near its poles. If the star spots 
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Table 5. Independent measurements of the mean density (p*), effective temperature (T cff ) and metallicity ([Fe/H]) for F1AT-P-11 and their weighted 
mean values. 


P*/Po 

Teff [K] 

[Fe/H] 

Source 

Notes 

i 7 c+0.92 
*' -0.39 

4780 ± 50 

0.31+0.05 

Bakos et al. (2010) 

Discovery paper 

2.42 + 0.10 



South worth (2011) 

Kepler Q0-Q2 

2 ' 13 -05 



Sanchis-Ojeda & Winn (2011) 

Kepler Q0-Q2 

2.55+0.11 

4792 ± 69 

0.33 + 0.07 

Deming et al. (2011) 

Torres et al. (2012) 

Kepler Q0-Q2 plus B-band and J-band 


4624 ± 225 

0.26 ± 0.08 

Mortier et al. (2013) 


2.52 ± 0.08 



Muller et al. (2013) 

Kepler Q0-Q6 


that are not occulted are equivalent to completely black spots 
that cover a fraction / of the stellar disc, then from equation (19) 
of Seager & Mallen-Ornelas (2003), we can estimate that the 

3 

density will be systematically too high by a factor (1 -/) *. The 
density of HAT-P-11 predicted by our stellar models based on the 
observed values of T e fr, T* and [Fe/H] assuming Ti so < lOGyr 
is p* = 1.9 ± 0.1p o . If this were to be explained by unocculted 
dark spots this would require / ~ 0.3. Such a large value of / 
can be ruled out in the case of HAT-P-11 because of the unusual 
orbit of the planet relative to the rotation axis of the star. The 
rotation and orbital axes are almost perpendicular and the ratio 
of the orbital and rotation periods is close to 1:6. This means 
the planet effectively scans the stellar disc along 6 different lines 
of longitude. Distortions to the light curve are seen due to the 
planet crossing star spots with a typical size of about 5° located 
in two bands of latitude, similar to the pattern of spots seen on 
the Sun (Sanchis-Ojeda & Winn 2011; Beky et al. 2014). These 
spots are too small to explain the large value of p* inferred from 
the transit light curve of HAT-P-11. 


In principle, we could find a set of stellar models that match 
the properties of HAT-P-11 by increasing the assumed mixing 
length parameter, but energy transport by convection is expected 
to be less efficient in magnetically active stars like HAT-P-11, 
not more efficient as a larger value of <+mlt would imply (Feiden 
& Chaboyer 2013). A more plausible explanation for the high 
density of HAT-P-11 is that this star has an abnormally high he¬ 
lium abundance. We used garstec to calculate a grid of stellar 
models identical to the grid of standard stellar models described 
above, but with the helium abundance increased by an amount 
AT = +0.05. With this “very helium enhanced” stellar model 
grid we find a good fit to the observed properties of HAT-P-11 
( x 1 - 0.34) for a mass (M*) = 0.74 ± 0.02M 0 and an age 
(t*) = 1.7 + 0.7 Gyr. Clearly, there is some additional system¬ 
atic error in these values because the actual helium abundance 
of HAT-P-11 is unknown. It is very likely that other stars show 
variations in their helium abundance around the simple linear re¬ 
lation between Y and Z that we have used. Another star in our 
sample with rj SOl b ~ 0 and y 2 ~ 0.5 for v = 0 degrees of free¬ 
dom is WASP-84. It is possible that this is a very young star for 
which the observational errors place the density slightly below 
the value predicted by our grid of standard stellar models, but 
the results for HAT-P-11 also opens up the possibility that this 
star has a more typical age and is helium-rich, particularly since 
it is not associated with any star forming region and does not 
show any other signs of extreme youth such as a high lithium 
abundance. 


4.2. Qatar-2 and inflated K-dwarfs 

It has long been known that some K-dwarfs appear to be larger 
by 5 per cent or more than the radius predicted by standard stel¬ 
lar models (Hoxie 1973; Popper 1997). This “radius anomaly” 
is correlated with the rotation rate of the star, but also shows 
some dependence on the mass and metallicity of the star (Lopez- 
Morales 2007; Spada et al. 2013). The dependence on rotation 
is thought to be the result of the increase in magnetic activity 
for rapidly rotating stars. Magnetic activity can affect the struc¬ 
ture of a star by producing a high coverage of starspots, which 
changes the boundary conditions at the surface of the star, or by 
reducing the efficiency of energy transport by convection. What¬ 
ever the cause of the radius anomaly in K-dwarfs, the existence 
of inflated K-dwarfs is one likely explanation for the appear¬ 
ance of stars in Table 3 that seem to be older than the Galactic 
disc (lOGyr, Cojocaru et al. 2014) or older than the Universe 
(13.75 Gyr, Hinshaw et al. 2013). 

One method that has been proposed to deal with the radius 
anomaly is to simulate the magnetic inhibition of convection 
by reducing the mixing length parameter (Chabrier et al. 2007). 
This phenomenological approach has some support from stellar 
models that incorporate magnetic fields in a self-consistent way 
(Feiden & Chaboyer 2013). Qatar-2 is one K-dwarf in our sam¬ 
ple that has an isochronal age that is clearly older than the age of 
the Galactic disc (P(t i so < lOGyr) = 0.001). We re-analysed the 
data for Qatar-2 using grids of stellar models similar to those de¬ 
scribed above but with various values of ctmlt from 1.22 to 2.32. 
We find that models with ctmlt ~ 1.4 can match the properties 
of this stars for ages less than lOGyr. 

The unknown value of o-vna is a source of systemic error that 
potentially affects the mass and age estimates for all the stars that 
we have studied. This includes the stars for which we have found 
a satisfatory fit to the observations for some age within the range 
expected (0 to 10 Gyr), i.e., it may be that we have assumed an 
inappropriate value of crMcr for some stars but that this is hidden 
by deriving incorrect (but plausible) values for the age and mass. 

4.3. Stars with consistent isochronal and gyrochronological 
ages. 

From Table 4 we see that there are nine stars for which there is 
good agreement between Ti s0 and r gyro (p T > 0.1). These include 
the stars HAT-P-11 and WASP-84 discussed above that may be 
helium-rich stars. As there is no good fit to the properties of these 
stars within our model grid we do not discuss them further in this 
section. 

Of the seven remaining stars, HD 209458 and CoRoT-4 both 
have reasonably precise gyrochronological and isochronal ages. 
For the other stars the agreement between the two age estimates 
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is partly due to the large uncertainty in one or both of these 
age estimates. The time scale for tidal spin-up of the star de¬ 
pends sensitively on the structure of the star and whether there 
is any resonance between the orbit of the planet and internal 
gravity wave modes in the star (Ogilvie 2014), so we should 
not expect that there is a precise value of logG^ai/yD that di¬ 
vides systems with and without tidal spin-up. Nevertheless, the 
good agreeement between the isochronal and gyrochronologi- 
cal age for HD 209458 and CoRoT-4 suggests that tidal spin- 
up for planet host stars may be inefficient for systems with 
log(Ttidai/yr) £ 12.5. 


4.4. Rapidly rotating stars 

Maxted et al. (2015) compared the observed masses of stars in 
detached eclipsing binaries to the masses predicted using bage- 
mass based on their density, effective temperature and metallic- 
ity. They found that masses of some stars with orbital periods 
less than about 6 days were under-predicted by about 0.15 M 0 . 
The rotation periods of these stars are expected to be equal to 
their orbital periods due to strong tidal interactions between the 
stars. There were no stars in the sample of eclipsing binaries used 
by Maxted et al. (2015) with orbital periods in the range 7-14 
days, so it is not known whether stars with rotation periods in 
this range are affected by the same problem. The cause of this 
discrepancy is not known so the isochronal ages for some stars 
with rotation periods less than about 6 days will be unreliable, 
and there is also the possibility that this problem affects stars 
with rotation periods up to about 14 days. 

There are four stars in our sample of planet host stars with 
rotation periods less than about 6 days (CoRoT-2, CoRoT-6, 
CoRoT-18 and Kepler-63). All four stars are clear examples 
of the gyrochronological age being significantly lower than the 
isochronal age ( p T < 0.02). Excluding these four stars from the 
sample does not have a strong affect on our conclusions as there 
are several other examples of stars with gyrochronological ages 
significantly lower than their isochronal age. 

If we take the more cautious approach and exclude all stars 
in our sample with rotation periods less than 14 days and HAT- 
P-ll (discussed above), we are left with a sample of 14 stars. 
Of these, about half have a gyrochronological ages significantly 
lower than their isochronal age (p T < 0.05). 

CoRoT-6 and Kepler-63 are examples of stars where the tidal 
forces between the star and the planet are very weak, e.g., weaker 
than the approximate limit log(r t id a i/yr) ^ 12.5 suggested in the 
previous section, above which limit it is reasonable to assume 
that there is no tidal spin-up for planet host stars. In these cases, 
we expect that the gyrochronological age of the star should be 
reliable. If this is the case, then it may be that these G-dwarfs 
stars are affected by magnetic inhibition of convection in a sim¬ 
ilar way to K-dwarfs. WASP-41 is another G-dwarf for which 
the gyrochronological age is significantly less than its isochronal 
age ( p T = 0.03) despite having a very long tidal time scale 
(log(r t idai/yr) = 12.3). Although this is not a rapidly rotating 
star (P rot = 18.4d) we mention it here because it is known to 
be a magnetically active G8-type star based on the appearance 
of chromospheric emission lines in its spectrum (Maxted et al. 
2011). This is consistent with the idea that some G-type stars 
have isochronal ages that are not reliable (too large) as a result 
of magnetic activity that is related to the rotation of the star but 
that is not directly caused by rapid rotation. 


4.5. Comparison to other studies 

Brown (2014) found a “slight tendency for isochrones to pro¬ 
duce older age estimates” but that the “evidence for any bias on 
a sample-wide level is inconclusive.” All 8 stars in that sample 
with directly measured rotation periods have been re-analysed 
here. For 7 of these stars, our results are consistent with those 
of Brown (2014). The exception is WASP-50, for which Brown 
uses a value of T e ff based on photometric colour (Brown, priv. 
comm.), rather than the lower and more accurate value that we 
have used here based on an analysis of the spectrum from Gillon 
et al. (2011). From a comparison of the gyrochronological and 
isochronal ages. Brown notes that of these 8 stars, 2 show “an 
age difference of a few Gyr”. We see very clear discrepancies 
between ri so and r gyro for 7 of these 8 stars, i.e., p T < 0.1, often 
much lower. We have restricted our analysis to stars with directly 
measured P mt values and increased the number of such stars 
studied to 28 so we very clearly see that ti s0 and T g y r0 disagree 
for the majority of stars common to both studies. A key differ¬ 
ence between our study and the study by Brown is the way that 
the uncertainties on the ages have been calculated. Brown cal¬ 
culated the range of isochronal ages corresponding to the “1-cr” 
error ellipse on T e ff and p* (the error on [Fe/H] was neglected) 
and then used this range as the standard deviation of a normal 
distribution to represent the posterior probability distribution for 
Ti S0 . The correlation between Ti so and r gyro via their mutual cor¬ 
relation with the assumed stellar mass was not considered. With 
our Markov chain method we can accurately account for the er¬ 
rors in Ti S o and r gyr0 from the errors in T e ff, p* and [Fe/H] and 
calculate the joint posterior distribution for these values (e.g., 
Fig. 3). This enables us to accurately calculate a statistic like p T 
that has greater statistical power than the method used by Brown. 

Pont (2009) used the full sample of 41 transiting exoplanet 
host stars known at that time to investigate whether there was 
empirical evidence for tidal evolution in these planetary systems. 
He identified two stars (HD 189733 and CoRoT-2) that show 
large excess rotation among the late-type stars in this sample 
(T e ff < 6000 K). This is consistent with the values of p T < 0.1 
that we have derived for both these stars. Additional constraints 
on the age of both these stars are now available, so we discuss 
these stars in more detail below. Pont’s study also includes the 
stars WASP-4 and WASP-5, but the only information on the ro¬ 
tation of these stars available at them time was rotation veloci¬ 
ties measured from spectral line broadening comparable to the 
instrumental resolution. With the rotation periods now available 
from photometry we find that both these stars show some evi¬ 
dence for excess rotation (p T < 0.1). CoRoT-4 is also present in 
both samples and shows no evidence for excess rotation in either 
study (p T = 0.4). 

4.6. Other age constraints 

Poppenhaeger & Wolk (2014) took a very different approach to 
testing whether hot Jupiters affect the rotation rate of their host 
stars. They compared the stellar coronal X-ray emission of 5 
planet host stars with their companion stars in wide binary sys¬ 
tems and used this activity indicator to estimate the age of both 
stars in each binary system. They found much higher magnetic 
activity levels in HD 189733 and CoRoT-2 than in their compan¬ 
ion stars. The estimated age of HD 189733 based on the X-ray 
flux is Tx = 1-2Gyr and for CoRoT-2 is rx = 0.1-0.3 Gyr, 
whereas their companions both have rx > 5 Gyr. In contrast, 
there was no evidence for a difference in age between the com¬ 
panion and the planet host star for t Boo (rx = 1 -2 Gyr), v And 
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(rx > 5 Gyr) or 55 Cnc (tx > 5 Gyr). There is very good agree¬ 
ment between our estimate of r gyro and rx for CoRoT-2 and T gyr0 
is consistent with rx in the case of 55 Cnc. The probability that 
Tiso for CoRoT-2 calculated using our grid of standard stellar 
models is consistent with the value of tx for its companion is 
P(rj so > 5 Gyr) = 0.09, i.e., if we assume that the value of rx for 
the companion is an accurate estimate of actual age of CoRoT-2 
then the isochronal age is likely to be an underestimate of the true 
age. Increasing the helium abundance by about +0.04 or reduc¬ 
ing <+mlt by about 0.4, or changing both these factors by about 
half as much, would be sufficient to bring the T; so and rx into 
agreement. For 55 Cnc, both Ti so and r gyr0 are consistent with 
each other and with the lower limit tx > 5 Gyr. Poppenhaeger 
& Wolk conclude from their observations that the presence of 
hot Jupiters may inhibit the spin-down of host stars with thick 
outer convective layers. With a sample of only 5 stars it is not 
yet possible to make precise estimates of what fraction of planet 
host stars may be affected by their planetary companions or what 
properties determine the strength of this interaction, but it is clear 
that r gyr0 is not a reliable estimate of the age for some planet host 
stars. 

The lithium abundance at the surface of a star can provide 
useful constraints on the age of a star, particularly for stars 
younger than about 600 Myr where there are good data in var¬ 
ious open clusters that can be used to calibrate the dependance 
of the lithium depletion rate with mass. We have compared the 
observed lithium abundances for the stars listed in the Table 4 to 
the calibration data from Sestito & Randich (2005) to see if the 
resulting constraints on the age are consistent with the values of 
Ti S o and r gyro . In general we can only set a lower limit to the age 
§ 0.6 Gyr that is consistent with both of these age estimates. A 
notable exception is CoRoT-2, which has a surface lithium abun¬ 
dance that implies an age for this star of ~ 0.2Gyr(Schroteretal. 
2011). This agrees very well with the gyrochronological age for 
this star. There does not appear to be any simple way to recon¬ 
cile the apparent young age of CoRoT-2 based on the available 
age indicators with the lack of X-ray emission from its K-dwarf 
companion. 

Another case where there is evidence to support the gy¬ 
rochronological age is WASP-77. Maxted et al. (2013) used 
the projected rotational velocity of the K-dwarf companion to 
WASP-77 to infer an gyrochronological age of 0.4' [j/ Gyr, con¬ 
sistent with their estimate of 1.0 /q 5 Gyr for the gyrochronologi¬ 
cal age of WASP-77. Maxted et ah (2013) used a different cali¬ 
bration for the gyrochronological age to the one used here, but it 
is clear that the gyrochronological age of companion to WASP- 
77 is inconsistent with the age for WASP-77 estimated using a 
grid of standard stellar models. If we assume a lower helium 
abundance for WASP-77 (Fig. 2) we will obtain a slightly lower 
estimate for the isochronal age, but it is not possible to reconcile 
the gyrochronological age and isochronal age without using a 
helium abundance less than the lower limit set by the primordial 
content of the Universe. It is possible to reconcile the isochronal 
age with the gyrochronological age by using a reduced value of 
the mixing-length parameter ctmlj much lower than the solar- 
calibrated value. WASP-77 shows chromospheric emission lines 
characteristic of magnetically active stars despite its moderate 
rotation period, P rot = 15.4d. Reducing the value of o-vilt is 
consistent with the idea of magnetic inhibition of convection that 
has been succesfully applied to explain the radius anomaly in K- 
dwarfs, but extends the idea in this case to a G-type star. 5 


5 We take the dividing line between K-type and G-type dwarfs to be 
0.8 M e in mass or T c fr = 5300 K. 


4.7. Reliability of gyrochronological ages for single stars. 

Meibom et al. (2015) have recently used Kepler photometry of 
30 stars in the open cluster NGC 6811 to measure their rota¬ 
tion periods. The calibration used here (Barnes 2010) predicts 
rotation periods for stars in good agreement with the observed 
values at the age of the cluster (2.5 Gyr) across the entire mass 
range studied (0.85-1.25 M 0 ). Epstein & Pinsonneault (2014) 
criticised the method used by Barnes (2010) to derive their cali¬ 
bration. We used Fig. 7 and Fig. 8 of their paper to estimate “by¬ 
eye” the gyrochronological ages of all the stars in our sample 
and found that these generally agree well with the values derived 
here, particularly if the calibration using the Kawaler style wind 
mass loss rate is used. For rapidly rotating stars the calibration 
of Epstein & Pinsonneault (2014) provides only an upper limit 
of 0.5 Gyr. This is consistent with the young ages that we derive 
for these stars, but it may be that these very young gyrochrono¬ 
logical ages are not reliable to the precision quoted here. This is 
unlikely to affect our conclusions substantially because an up¬ 
per limit of 0.5 Gyr is generally sufficient to show that these 
rapidly rotating stars have gyrochronological ages inconsistent 
with their isochronal ages. The calibration of Epstein & Pinson¬ 
neault (2014) was also applied to the data for stars in NGC 6811 
by Meibom et al. (2015) and found to give rotation periods that 
are too low by about 10% for stars with masses 0.85 M 0 but 
that agree well with the observed rotation periods for stars with 
masses near 1 M 0 . 


5. Conclusions 

By using new data and improved analysis methods we have 
shown that there is now good evidence that some exoplanet 
host stars rotate more rapidly than expected. For our sam¬ 
ple of 28 transiting exoplanet host stars, about half the sam¬ 
ple have gyrochronological ages that are significantly less than 
their isochronal ages. In a few such cases there are indepen¬ 
dent constraints on the age of the star that are consistent with 
the isochronal age, which suggests that tidal spin-up of the host 
star has occured in these systems. However, in several cases it is 
not clear that tidal interactions between the star and the planet 
are responsible for this discrepancy. For some K-type stars this 
is a result of the well-known radius anomaly that may be due to 
magnetic inhibition of convection. We find some evidence that 
this anomaly may also affect some of the rapidly rotating and/or 
magnetically active G-type stars in our sample, either from in¬ 
dependent age constraints on the age of the star (WASP-77) or 
because it may be implausible that the strength of the tidal in¬ 
teraction with the planet is strong enough to spin-up the star 
(CoRoT-6, Kepler-63 and, perhaps, WASP-41). 

Some planet-host stars (HAT-P-11 in particular) appear to be 
much denser than predicted by stellar models. These stars may 
be significantly enhanced in helium. This makes it difficult to 
assess the reliability of isochronal mass and age estimates for 
these stars. There is currently no simple explanation for the in¬ 
consistency between the young age of CoRoT-2 implied by stel¬ 
lar models, gyrochronology, its X-ray flux and its high lithium 
abundance with the very old age inferred for its K-type compan¬ 
ion based on its lack of X-ray flux. 

Our improved analysis methods have enabled us to show that 
there is now clear evidence that the gyrochronological ages of 
some transiting exoplanet host stars are significantly less than 
their isochronal ages. However, a careful consideration of all the 
available data on a case-by-case basis shows that it is not always 
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clear that this is good evidence for tidal spin-up of the host star 
by the planet. 
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